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Abstract. As stars close to the galactic centre have short orbital periods it has been possible to trace large fractions of their 
orbits in the recent years. Previously the data of the orbit of the star S2 have been fitted with Keplerian orbits corresponding to 
a massive black hole (MBH) with a mass of Ms/f=3-4xlO*Mo implying an insignificant cusp mass. However, it has also been 
shown that the central black hole resides in a ~r' diameter stellar cluster of a priori unknown mass. In a spherical potential 
which is neither Keplerian nor harmonic, orbits will precess resulting in inclined rosetta shaped trajectories on the sky. In 
this case, the assumption of non-Keplerian orbits is a more physical approach. It is also the only approach through which 
cusp mass information can be obtained via stellar dynamics of the cusp members. This paper presents the first exemplary 
modelling efforts in this direction. Using positional and radial data of star S2, we find that there could exist an unobserved 
extended mass component of several IO'Mq forming a so-called 'cusp' centered on the black hole position. Considering 
only the fraction of the cusp mass Mi-2„,,„ within the apo-center of the S2 orbit we find as an upper limit that Ms2„^„/(Mbh + 
Ms2„p„) < 0.05. A large extended cusp mass, if present, is unlikely to be composed of sub-solar mass constituents, but could 
be explained rather well by a cluster of high M/L stellar remnants, which we find to form a stable configuration. 
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1. Introduction 

Over the last decade, evidence has been found for the exis- 
tence of massive black holes (MBH) in the centres of many 
nearby galaxies. With increasing observational data - 12 
MBHs candidates detected until 1995 (Kormendy & Rich- 
stone 1995), and more than 37 until 2001 (Kormendy 2001, 
Ferrarese et al. 2001) - it is argued that most galaxies harbor 
nuclei dominated by MBHs with masses that range between 
10^ and 10^^ solar masses. 

Located at a distance of only ~8 kpc from the solar sys- 
tem (Reid 1993; Eisenhauer 2003), the Galactic Centre (GC) 
is the closest and therefore best object for investigating physi- 
cal processes in the galactic nucleus of a typical spiral galaxy. 
It offers a unique "laboratory" for studying stars and gas 
in the sphere of influence on a super-massive black hole, 
(e.g., Genzel, Hollenbach, & Townes 1994; Morris & Ser- 
abyn 1996; Mezger, Duschl, & Zylka 1996; Melia & Fal- 
cke 2001), with a degree of detail that cannot be accessed 
in any other galactic nucleus in the foreseeable future. With 
the now available high sensitivity and angular resolution. 
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large ground-based telescopes offer the opportunity to obtain 
an unprecedented view of the Galactic centre. Initially, with 
speckle imaging techniques and lately with adaptive optics 
techniques, high angular resolution images on the Galaxy's 
central cluster were obtained. This first set of observations 
was able to measure stellar motions on the plane of the sky, 
yielding estimates of the projected velocities (Eckart & Gen- 
zel 1996; Ghez et al. 1998), projected accelerations (Ghez et 
al. 2000; Eckart et al. 2002a), and three-dimensional orbital 
motions (Eckart et al. 2002a, Schodel et al. 2002, 2003; Ghez 
et al. 2003b), which each provided a successively stronger 
case for a super-massive black hole at the centre of the Milky 
Way and its association with the unusual radio source Sgr A* 
(Lo et al. 1985). 

More than a decade of high-resolution infrared observa- 
tions of proper motions in the GC, with the ESO New Tech- 
nology Telescope (NTT) and the ESO Very Large Telescope 
(VLT) (Eckart & Genzel 1996; Eckart et al. 2DD2a Schodel 
et al. 2002, 2003), as well as with the Keck telescope (Ghez 
et al. 1998, 2000, 2003b), have revealed at least 6 stars that 
show substantial acceleration due to the super-massive black 
hole associated with Sgr A* and are on bound orbits around 
it. With a peri-centre of less than 0.6 mpc (15 mas) and an 
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orbital period of ~ 15 years, S2 is the most striking case of 
these. 

A series of observations with the NAOS/CONICA adap- 
tive optics system/near-infrared camera at the ESO VLT unit 
telescope 4 that covered the peri-centre passage of the star 
S2 around Sgr A* allowed Schodel et al. (2002) to approxi- 
mate a Keplerian orbit and to measure the enclosed dark mass 
down to a distance of ~ 0.6 mpc from Sgr A*. With these ob- 
servations, they could exclude a neutrino ball scenario (Mun- 
yaneza & Viollier, 2002) as an alternative explanation for the 
dark mass concentration. They excluded as well a cluster of 
dark astrophysical objects (Maoz 1998) such as neutron stars, 
leaving a central super-massive black hole as the most prob- 
able explanation. 

Using the Keck 10m telescope, Ghez et al. (2003a) mea- 
sured a Keplerian orbit for the star S2 which agreed with 
the results of Schodel et al. (2002). Additionally, Ghez et 
al. (2003a) reported the first detection of spectral absorption 
lines (both Bry (2.1661 fim) and He I (2.1126 jum)), provid- 
ing the first line-of-sight velocity measurements of star S2. 
These measurements resolved the ambiguity on the inclina- 
tion of the S2 orbit indicating that its position was behind the 
black hole when it passed through its peri-centre. In addition, 
stellar rotational velocities suggest that S2 is an O8-B0 dwarf 
star and thus a massive (~ 15 M©) young star (< 10 Myrs). 

From data taken with NAOS/CONICA and the new NIR 
integral field spectrometer SPIFFI at the ESO VLT, Eisen- 
hauer et al. (2003) reported new astrometric observations and 
additional spectroscopic observations of the star S2, reducing 
the uncertainties on the orbital parameters. They also gave the 
most accurate primary distance measurement to the centre of 
the Milky Way of 7.94 ± 0.42 kpc, which is in agreement 
with earUer determinations (see Reid 1993). 

In this paper, we explore the possibility that there exists 
a compact, continuous mass distribution composed of several 
undetectable faint stars or perhaps some more exotic mate- 
rial in addition to the point mass of Sgr A*. In this case, the 
orbit of S2 will not follow a simple Keplerian orbit, but will 
rather show peri-centre-shifts that result in rosetta shaped or- 
bits. This idea is motivated by the observation that the steUar 
density does not flatten out, but exhibits a steep peak towards 
the centre, a so-called cusp (Eckart et. al, 1995; Alexander , 
1999; Genzel et al., 2003). In contrast to earlier studies, the 
main approach of this work is that the mass-to-light ratio, 
M/L, is not considered to be constant over the entire range 
of the GC stellar cluster. Indeed, there are indications that the 
stellar population varies with position and is not quite well 
mixed (Alexander, 1999). Furthermore, the exact composi- 
tion of the cusp is still unknown and our current understand- 
ing of the stellar distribution in the GC is incomplete (Genzel 
et al. 2003). With current observations, low-mass stars iK> 
21 mag) cannot be observed in these dense cusp regions and 
the true value of the M/L ratio is not known. Also, as pointed 
out by Baumgardt et al. (2003), dynamical evolution of a 
dense stellar cluster will result in a strong increase of M/L ra- 
tio by segregation of stellar evolution remnants to the centre. 
On the other hand, the approximation that the dynamics in the 
central region is Keplerian is directly related to the implicit 



assumption of Eckart et al. (2002a), Schodel et al. (2002), and 
Ghez et al. (2003a) that the M/L ratio at 2 /urn in the cusp is 
as low as in the outer stellar cluster (M/L = 2Mq/Lq). 
Considering these circumstances, stars with short orbital pe- 
riods, in particular the star S2, of which orbital data are best 
known, play a key role in exploring the gravitational poten- 
tial. We show that the present observational data on S2 cannot 
discriminate between a Keplerian and a non-Keplerian po- 
tential. Subsequently, we study the influence of an extended 
distribution of dark mass near Sgr A* taking into considera- 
tion the constraints set by the orbit of S2 as well as the limits 
set on the total enclosed mass at larger radii. In this study 
relativistic effects are neglected as they are second order cor- 
rections. 

In § 2 we outUne our modelling of the density distribu- 
tion of the central region and show that the star counts near 
Sgr A* can be described by a compact Plummer model core. 
The method used to compute orbits in non-Keplerian poten- 
tials is explained in§ 3. In § 4 possible orbital models for S2 
are presented. We discuss the implications of the limits on the 
cusp mass derived from the non-Keplerian orbit modelling in 
§ 5, and draw our conclusions in § 6. 

2. Modelling the stellar cluster and the enclosed 

mass 

In this section it is described how the most recent observa- 
tions of the GC star cluster (Genzel et al., 2000, 2003) can be 
used to model the mass distribution function as close as a few 
milli-parsecs from Sgr A*. A surface density profile is fitted 
to number density counts and used, together with dynamical 
mass estimators, to deduce the 3 dimensional mass distribu- 
tion function. By following this approach, we have assumed 
that the visible stars trace the mass distribution of the cluster. 

2.1. Density model of the stellar duster 

Genzel et al. (2003, see their Fig.7) determined the surface 
number density of stars in the central parsec by counting 
the sources in annuU with increasing radius around Sgr A*. 
Within ~ 2 arc-seconds of Sgr A*, their data indicate an ex- 
cess of faint stars above the counts expected in a cluster with 
a flat core, as it had been already suggested earlier (Eckart et 
al. 1995; Alexander 1999). They confirmed the existence of 
a "cusp" and showed that it is centred at a position of (A a, 
K6) = (+0.09", - 0. 15") from Sgr A*, wifli an uncertainty of ± 
0.2". Here we use the same K-band counts {K< 18) together 
with the earlier SHARP/NTT counts {K < 15), without taking 
into our study their higher magnitude Umit H-band counts {H 
< 19). In contrast to the 'by-eye' H-band counts, the K-band 
counts, were determined by an automated procedure, they are 
therefore possibly less biased, and cover the whole range of 
the cluster, not only the central arc-second. 

Genzel et al. (2003) fitted a broken power-law to their 
stellar counts, with a distinct break at 10 arcsec. In this sec- 
tion, we show that it is also possible to fit the number density 
counts with a superposition of several Plunnmer models. The 
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Fig. 1. Surface density of stars vs. distance from Sgr A*. 
The downward pointing triangles represent the CON- 
ICA/NAOS data for Ks<18, the upward pointing triangles 
are SHARP/NTT data for Ks<15 stars scaled up by a fac- 
tor of 5 in order to match the fainter NAOS/CONICA counts. 
The data are taken from Genzel et al. (2000, 2003). The long- 
dashed line is the superposition of the 5 outermost Plummer 
models (see. Table The addition of the innermost Plum- 
mer model (the 'inner' cusp') gives the 'composite model'. It 
is represented by the 3 other curves, where the thick-straight 
line corresponds to the middle values of the 'inner cusp' pa- 
rameters Ri and Mi, the thin-straight line to its lower values, 
and the short-dashed line to its upper values, and the short- 
dashed line to its upper values. 

Plummer model has the advantage of being analytically in- 
tegrable. It resembles actual clusters with compact cores and 
an extended outer envelope (Spitzer 1975). 

In this model the 3 dimensional radial density distribution 
p{r) and the projected mass density cr{r) are described by: 

' (1) 



P(r) = 



4^^^1 + ^)^ 



4-p(Q)R 



1 



(2) 



3 (1 + ^)2 
where R is the core radius and M the total mass. 

Our best surface density fit to the data consists of a super- 
position of 6 different Plummer models, which will be termed 
in the following as 'composite model': 

' M, R} 



cr(r) 



(3) 



,=1 - (r'-+Rjr 

The six values of the core radii, are given in Tabled 
Fits with larger numbers of Plummer models are possible as 



Table 1. The parameters of the 'composite model' used to 
fit the GC stellar cluster. For each Plummer model we give 
the value of the core radius R, (in arcsec and in parsec) de- 
duced from the surface density plot. M, is the corresponding 
total mass scaled to the enclosed mass distribution (Fig. |2jl. 
Mgpoi is the contribution in mass inside a sphere which has a 
radius equal to the apo-astron of S2. p, is the calculated aver- 
age density in Mg.pc"^. The first Plummer model (Ri = 0.4" 
or 0.015 pc) corresponds to the inner cusp centered at Sgr A*. 
The addition of the remaining 5 models represents the density 
distribution of the outer stellar cluster. 
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well; six was the smallest numbers of components with which 
we were able to represent the data, i.e. the cusp plus the stellar 
cluster. 



2.2. The enclosed mass 

Similarly to eq|3 the mass distribution of the 'composite 
model' is given by 



M(r) = 



^(,2+^2)3/2 



(4) 



Under the assumption that the visible stars trace the mass, 
the M, can be determined by fitting eql^to the dynamical 
enclosed mass estimates. 

In order to calibrate the enclosed mass M(r) we only rely 
on a solid measure of the central mass (within 0.001 pc from 
SgrA*) derived from the Keplerian orbit of the star S2 and 
mass estimates of the entire stellar cluster at distances of 
several parsecs from the center At these large radii the pub- 
lished results of all different statistical mass estimators give 
very consistent results within uncertainties that are fully ac- 
ceptable for our study (Eisenhauer 2003, Genzel et al. 2000, 
Serabyn & Lacy 1985; Gusten et al. 1987; Lacy et al. 1991; 
Herbst et al. 1993; Roberts & Goss 1993). By tying our mass 
modelling as a function of radius to the stellar number den- 
sity distribution, it does not depend on possible anisotropics 
of the velocity field at radii of less than a parsec. It is only 
there where the anisotropics are apparently present at a con- 
siderable level (Genzel et al. 2000, Schodel et al. 2003). In 
particular we would like to point out that we are not fitting 
to enclosed mass estimates at these small radii. In Fig. 2 these 
are only shown for completeness. The only assumption that 
goes into our considerations is that of a constant M/L ratio 
as a function of distance. This assumption is relaxed and dis- 
cussed in the following sections of the paper. 

Our best fit to these dynamical mass estimates corre- 
sponds to a 2.9x10*' Mo + 0.3x10'' Mq black hole plus the 
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mass distribution of eq|4l of which the 6 total masses M,, are 
hsted together with their corresponding densities, p,, in Ta- 
ble[2 An offset in the value of the enclosed compact mass will 
not influence our modelling significantly because the shape 
and mass of our cluster composite model depends primarily 
on the profile of the stellar number counts. 

The combination of the five outermost components gives 
a similar configuration as the flattened isothermal sphere of 
core radius ~ 0.3 pc derived by Genzel et al. (2003) (see 
mass densities below). Fig.^shows that there exists a mini- 
core centered on Sgr A*, which we denominate by the 'inner 
cusp' . It is represented by the innermost Plummer model with 
a core radius Ri = 15.5mpc = 0.4". 

With an apo-centre of approximately 9 mpc (Schodel, 
2004; Ghez et al., 2003a), the star S2 orbits in a region in- 
side the core radius (7?i~15mpc) of this innermost model. 
Therefore its path can be mainly influenced, in addition to 
the BH potential, by the 'inner cusp' gravitational potential 
(see colomn 5 in Table0. At these distances from the centre, 
the present number density counts have large error bars (see 
Fig.Q. The limits on the core radius, Ri, of our model, de- 
duced from the error bars on the number density counts, are 
found to be 13.2 mpc for the lower value and 20.2 mpc 
for the upper one R""'\ Varying the mass. Mi, of the central 
Plummer model component allows to study cases in which 
the mass-to-light ratio, M/L, varies as a function of distance 
from Sgr A*. As described in § 3 and § 4, exploring the pos- 
sible ranges foiRi and Mi that fulfill the constraints given by 
the orbital measurements of star S2 allows to derive dynami- 
cal information on the cusp mass. In the following, we make 
use of this model which fixes the density profile of the cusp 
to the stellar profile. 

In Fig.|2lthe long-dashed line corresponds to our best fit, 
the thick-straight line shows only the enclosed mass without 
a BH obtained with our model of the extended cluster The 
thin-straight and the short-dashed lines represent the above 
discussed extreme cases of an inner cusp radius Ri - R""" = 
13.2 mpc with M = 85400 and Ri = R"^"''' = 20.2 mpc with 
M - 8820q. Note that these errors result from the uncertain- 
ties of the stellar number counts at the smallest distances to 
Sgr A* (Fig.dl- 

The density of the 'inner cusp' is a few hundred times 
larger than that of the second Plummer model and it has a 
~ 33 times smaller core radius (see Fig. Hand TableQ. Our 
model implies density values as high as 1.68x10** Mopc^^* at 
0.1", and 1.54x10^ Mopc^^* at 1". These values are compara- 
ble to the ones derived by Genzel at al. (2003) from a broken 
power-law density profile (7x10** Mopc"^ at O.T'and 3x10^ 
Mopc"^ at 1"). 

In the following we will use the above derived mass dis- 
tribution to determine the path of S2 in response to this non- 
point mass potential plus the central black hole. 

3. Modelling the orbit of S2 

In order to fit non-Keplerian orbits to the measured time- 
dependent positions of the star S2 in the above described type 
of mass potential, one needs a high accuracy integrator For 
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Fig. 2. Mass distribution in the central 10 pc of the Galactic 
Centre for an 8 kpc distance (Reid et al. 2003; Eisenhauer 
et al., 2003). The star point denotes the mass derived from 
the orbit of S2. The error bar combines the orbital fit and 
astrometry errors (Eisenhauer et al. 2003). Filled circles are 
mass estimates from a parameterised Jeans-equation model 
of the radial velocities of late type stars, assuming isotropy 
(Genzel et al. 2000, Schodel et al. 2002). We omitted ail 
mass estimators at radii of less then 1" since they are effected 
by the surface mass density of the cusp which had not been 
taken into account in previous investigations. Filled upward- 
pointing triangles denote Leonard-Merritt projected mass es- 
timators from a new NTT proper motion data set (Schodel 
et al. 2002). Filled down-pointing triangles denote mass es- 
timates from the ionized and neutral gas dynamics (Serabyn 
& Lacy 1985; Gusten et al. 1987; Lacy et al. 1991; Herbst 
et al. 1993; Roberts & Goss 1993). A comparaison of repre- 
sentative mass models to the measured velocity dispersions is 
given in Genzel et al. (1996, 2000). The long-dashed line is 
a fit to the measured enclosed mass against distance. That re- 
sults from a combination of the mass distributions (see § 2.2) 
deduced from the composite model fit to the stellar cluster, 
plus a 2.9x10^ M0 black hole. The same model without a 
point mass is represented by the 3 other curves, they corre- 
spond to acceptable models of the stellar cluster within the 
error bars. 

this purpose, we chose a fourth order Hermite integrator de- 
rived from the one used in high-accuracy A^-body simulations 
(Aarseth 1999, Spurzem 1999, for an introduction of the Her- 
mite scheme see Makino & Aarseth 1992). The advantage 
of the Hermite scheme is that it allows a fourth order accu- 
rate integration based on only two time steps. Therefore it re- 
quires the analytic computation of the time derivative of the 
gravitational force. This is the point where the integrability 
of the Plummer model is convenient. We computed the tra- 
jectory of S2 around the BH and through the extended mass 
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Table 2. Characteristic ;i'^tof of our fits. The first column lists 
the obtained characteristic for which the resulting orbits 
are considered to be fitting orbits. The 2'"' column expresses 
the total non-reduced ;if^fof for 33 degrees of freedom (x^roi = 
X^totP^)- When column 3 lists the deviation of the obtained 
xjo, from the lowest obtained column 4 corresponds their 
resulting cr error bars. 
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as given by our mass model (see § 2). The presence of an ex- 
tended mass induces retrograde pericentre-shifts which result 
in open rosetta shaped orbits. 

In total there are ten fit parameters - six orbital elements, 
the total central mass M,,,,, the fraction / of the extended mass 
to the total mass M,„,, and the 2-D position P of the central 
mass. All these values are known within a certain margin. 
We used the position measurements for S2, the dynamical 
position of Sgr A* and the 5 line-of-sight velocity measure- 
ments as determined by Schodel et al. (2003), Eisenhauer et 
al. (2003) and Ghez et al. (2003a). 

A valid computed orbit has to fit both the measured ve- 
locities and positions. Therefore, two reduced values, one 
for the positions, x^pos, and another for the velocities, x^vei, 
were estimated. While fitting the orbits, we considered the 
X^pos values for the positions and velocities independently 
because the x^ values for the velocities were, in general, rel- 
atively much smaller than for the positions. In some cases 
they indicated that measurement errors were systematically 
overestimated. This way, it was possible to attain, both in 
velocities and in positions, a similar fit quality. Under this 
condition, it is always possible to define a common by 
adding both the x^ values obtained from the positions and the 
velocities. The total reduced x^, xfot^ simply the x^tot di- 
vided by the number of degrees of freedom which is equal to 
33 in our case. Table |2l lists obtained total reduced and 
non-ieduced x^ tot which indicated errors < 5cr. These errors 
were estimated by applying a deviation Ax^tot from the low- 
est obtained ;^'^,o,. In our following analysis we used orbital 
fits with an error < Icr. 

3.1. Influence of model parameters 

In order to study the main effects of a varying cusp mass 
on the orbit of the star S2 we investigated a representative 
coverage of the parameter space that includes the full range 
of possible scenarios. In each single orbital run, the total 
central mass, M,o, and its position as well as the extended 
mass fraction, /, were fixed; only the 6 orbital elements 
are varied. The three parameters, M,o,, the central mass 
position and /, were varied independently forming a fully 
representative 3-D grid. In total ~360 orbits were computed 



and fitted to the data. In the following, we give a short 
description of this parameter space: 

The total mass M^/, together with the the three veloc- 
ities and three coordinates determine the orbital elements. 
From the imaging data, sky positions and proper motions are 
given in a certain error range for each epoch. For a chosen 
initial epoch, these six input values are varied within their ob- 
servational uncertainties. Our minimization method is a grid 
scheme, which, although expensive in computing time (time 
oc N^, N being the number of steps per parameter), is a robust 
method which does not fall easily on local minima. Choosing 
an appropriate for each parameter, x^ values were com- 
puted and a 6-D grid was constructed. Ultimately, the lowest 
X^ value determines the best fitting orbit. 

The black hole is treated as a point mass with the ex- 
tended component being centered on it. The lower central 
mass limit (Reid et al. 2003) and upper central mass limit 
(Ghez et al. 2003; Schodel et al. 2003) are well determined. 
Based on these earlier works, the value of the total central 
mass is varied between 2.6 xIO^'Mq and 4.8 xlO''Mo, consid- 
ering six different values of M,o,. In addition, for any of these 
given total central masses values, M,,,,, the effect of various 
extended mass fractions, /, was investigated (see Fig.^. 

We used the position of Sgr A* as determined by Reid 
et al. (2003) with an uncertainty of 10 mas. The location of 
the central mass was varied at 14 different positions within 
this error bar. Fig. |3l shows these positions along with their 
obtained x^tot values. These were fairly distributed around 
the offset position which is the deviation of the focus of the 
Keplerian orbit from the central radio position the way deter- 
mined by Eisenhauer et al. (2003). In order to check the accu- 
racy of the focus of the Keplerian orbit, we compare the two 
positions for the star S2 at the epoch 2002.30 calculated with 
the orbital parameters given by Eisenhauer et al. (2003) and 
by Ghez et al. (2003b). From the corresponding positions and 
uncertainties, an offset of the projected orbital focus given by 
Ghez et al. (2003) from the one by Eisenhauer et al. (2003) 
of 0.0 ± 2.9 mas East and 0.8 + 2.0 mas North can be derived. 
Hence, we find an excellent agreement between the two inde- 
pendent results. 

In order to take into account the uncertainty of the core 
radius of the inner cusp, the orbital fits computed for our 
best fitting model value of Ri - 15.5 mpc, were repeated 
only in the case of the fitting orbits with an error < 2 cr, for 
the extreme cases of R"^'" = 13.2 mpc and R'""^ = 20.2 mpc. 

4. Results of non-Keplerian fits 
4.1. Orbits with lowest ^^,o, values 

Our results show that the available data on the path of S2 can 
be fitted very well with non-Keplerian orbits. However, at 
present the orbit of S2 is not sufficiently well determined to 
exclude a pure Keplerian orbit (see Fig. |6j. In fact, in both 
cases a minimum ^^,o, value of ~0.7 was obtained. In other 
words, it cannot be excluded that there exists an extended 
mass component in addition to the black hole. In the case of 
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4.1 Orbits with lowest x^ioi values 
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Fig. 3. The different positions of the central mass as they were 
considered associated with the lowest obtained tot value at 
each given position. The offset position is centred in the area 
showing the lowest Jp^fo, values. 

Keplerian orbits and for different central mass position, fits 
for a 3.65xlO^M0 black hole shows the lowest x^,o, value. 
If we assume an additional extended mass component, the 
lowest x^tot values corresponded to fits for a total mass of 
3.65xlO^M0 of which 5% are extended or a 4.1xlO''M0total 
mass of which 10% or 15% are extended, or 4.45xlO^M0 to- 
tal mass of which 20% are extended, as well as a 4.8xlO^M0 
total mass of which 25% are extended. We can also consider 
only the fraction of the cusp mass Ms2„^„ which is located at 
distances from the MBH smaller than the apocenter of the S2 
orbit. In this case the 25% limit corresponds to 

Ms2„,J(Mbh + Ms2„,J< 0.05. 

The orbital parameters for these < 1 cr fitting orbits are given 
in Appendix A together with the full parameter range that 
produces orbits with an error of <5cr . Within the uncertain- 
ties the orbital parameters for the Keplerian case agree with 
earlier work (Schodel et al. (2003), Ghez et al. (2003), Eisen- 
hauer et al. (2003)). 

4.1.1. Position of Sgr A* 

Only orbits with focuses at a, b, c, d and offset positions (see 
Fig. Hand Fig.|5jl resulted with an error of <lcr. This gives 
a constraint on the SMBH black hole position of only a sep- 
aration of +0.05 mpc (+1.3 mas) in R.A or in Dec. from the 
offset positions. 

4.1.2. Higher total central masses and cusp masses 

Fig.lHshows a plot of x^ tot vs. different discrete central total 
mass values and for each fraction / of the extended mass at 
the offset position. We also examined orbits with total masses 
higher than 4.8x10*' M0, up to 1 LexlO*" M0. We found that 
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Fig. 4. Plot of the x%t values vs. the total central mass, M,o,, 
calculated for different extended mass fractions, /, with the 
point mass at the offset position. Values of the same / are 
fitted with splines. Crosses: / = 0; "x" symbols: / = 0.05; 
stars: / = 0.1; empty squares: / - 0.15; filled squares: 0.2, 
empty circles: / = 0.25; and filled circles: / - 0.3. The two 
horizontal lines denoted by <1 cr and 5 cr delimit the region 
of fitting orbits (x^^, < 1.44). 
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Fig. 5. The lowest obtained reduced^^^; values in function of 
the separation from the offset position. These are fitted with a 
function of the form 7.2x^ + 0.5x + 0.7. 



at a 5 cr level of confidence the S2 orbit could fit the data 
for cusp masses of up to 70% of the total mass. This resulted 
in an upper limit of cusp masses as high as 8x10* cor- 
responding to an extended component of ~lxlO* M0 inside 
the sphere of a radius equal to the apocentre of S2. However, 
the total masses in these cases are in disagreement with the 
dynamical mass estimates at large radii. Therefore, we regard 
them as unrealistic but mention the result for completeness. 

Large amounts of extended mass inside the S2 orbit 
would be possible if the mass distribution function at the cen- 
tre were steeper than in our model or if the core radius were 
smaller than the apo-centre of S2 (~ 9 mpc). In that case we 
found that mass densities - exceeding the currently derived 
few 10** Mopc"^ by three or four orders of magnitude - are 
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Fig. 6. Three exemplary orbits determined by our fitting. Upper left panel: Keplerian orbit with 3.65xlO^M0 point mass. Up- 
per middle panel, non- Keplerian orbits with 4.1xlO^Mo total mass, thereof 10% extended. Upper right panel, non-Keplerian 
orbits with 4.8xlO^Mo total mass, thereof 25% extended. Here the central mass is at the offset position, 0.09mpc east and 
0.12mpc south from the nominal radio position of Sgr A*. The lower panels show the velocity in the relative R.A, the relative 
Dec, and along the line-of-sight as a function of time for the case of 4.1xlO^M0 total central mass with a 10% extended 
component. This case is representative for the other cases. The direction of A a the pericentre-shift is shown by an arrow. 



unlikely to form a stable configuration because stellar colli- 
sions will become of increasing importance (see section 5.2). 

4.1.3. Core radius and slope of the density profile 

In the case of the offset position, for all the orbits with errors 
< 2 cr, the fitting procedure was repeated for the extreme 
cases of /?™"=13.2 mpc and /?'["'^=20.2 mpc of the inner 
cusp radius. Even with these different core radii, the case of 
the offset position is found to give the lowest x^iot values. 
Also the fit quality remains unchanged at x^,o, - 0.7. We 
also consider that the compact component can have a much 
smaller core radius than the stellar profile, possibly well 
within the S2 apocentre. In this case, following the result 
of Rubilar & Eckart (2001, see their Fig. 5), we deduce 
that for the same extended mass fraction, /, the induced 
pericentre-shift is smaller for a bigger core radius. Thus for 
a core radius that does not exceed a value 15.5 mpc, our 
estimation on the mass of the cusp is an upper limit. 

Using a series of Plummer model density profiles, 
implies a value of the cusp slope of o- < 5. As mentioned in 
§ 2, Genzel et al. (2003) fitted the stellar number density with 
a broken power-law of a slope equal 1.4±1. From theoretical 
models of cusps around massive black holes, one would 



expect power-law indecies in the range 3/2 < 7/4 (BahcaU 
& Wolf 1977; Murphy et al. 1991). To study the influence 
of the denstiy profile on the pericentre-shift, we performed 
fittings of the S2 orbit considering a power-law profile of 
a slope equal to 1.4. The pericentre-shift induced by the 
same amount of mass inside S2's apocentre but different 
fractions of mass, /, did not exceed the one induced by a 
Plummer model density profile by more than 12%. Thus, by 
covering a wide range of slopes 1.4 <a< 5, we found that 
for the shallower mass distribution, the pericentre-shift never 
exceeds 12% of our initially derived value, thus our estimates 
indicate an upper limit on the mass that could be present 
inside the S2 apocentre, and consequently on the mass of the 
cusp. 



5. Discussion 

In the previous sections it was demonstrated that the mea- 
sured positions and line-of-sight velocities of S2 are not sufii- 
ciently well determined to exclude a Keplerian orbit. Assum- 
ing a simple Keplerian case, the central gravitational potential 
by definition is entirely dominated by the point mass associ- 
ated with Sgr A* and no dynamical constraints can be derived 
on any extended mass component due to e.g. the surrounding 
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5. 1 Extending the K-band Luminosity Function to Faint Stars 



stellar cluster. More revealing, however, is the more physical 
assumption of non-Keplerian orbits. The presence of a stellar 
cusp shows that - at least to a certain extent - there is some 
extended mass present near the black hole. This potential can 
be modelled by a central black hole plus an extended mass 
component. We assumed that the cusp has the same shape as 
the observed stellar number density counts which we fitted 
by a Plummer density model (see §2). We also took into ac- 
count that the total enclosed mass at larger radii (Ipc) would 
not exceed a value of ~ 4.8 x 10^ Mq (see section 3.1). 

Having shown that the current data allow for an extended 
mass component higher than the one of 8540 inferred 
from stellar number density counts which were scaled to 
match the velocity distribution at large radii (see §2), we 
will discuss its nature in sections §5.1 and §5.2. The value 
of 8540 Mq, representing the total mass of the 'inner cusp', 
is equivalent to a mass of 3100 Mq inside the core radius Ri 
(15.5 mpc or 0.4"). It is based on the assumption that the ra- 
tio of the stellar number counts to total stellar mass does not 
vary with radius and environment in the Galactic Centre, i.e. 
that the mass-to-light ratio M/L(2//m) is constant (see also 
Genzel et al., 2003). It is very likely that this assumption is 
unjustified. Genzel et al. (2003) find that the stellar popula- 
tion in the cusp differs to a certain degree from the popula- 
tion of the surrounding, large-scale cluster. Also, effects such 
as mass segregation and stellar collisions might work very 
effectively in the dense environment of the cusp and M/L 
would show higher values towards the centre (Baumgardt et 
al. 2003). Therefore, we consider that the mass-to-light ra- 
tio, M/L (2jum) of the 'inner cusp', can be different from that 
of the outer stellar cluster In the following, two hypotheses 
will be discussed: The existence of a cluster of faint low-mass 
stars, not yet detectable with the resolution and sensitivity of 
current instruments, and the existence of a cluster of heavy 
dark objects like stellar black holes and neutron stars. 

5.1. Extending the K-band Luminosity Function to Faint 
Stars 

We use the /iT-band luminosity function (KLF) for K< 18 mag 
within a projected radius of 1.5" from Sgr A* as it was de- 
termined by Genzel et al. (2003 - see their Fig.l 1). Note that 
the KLF is not corrected for extinction or for the distance 
modulus. The authors fitted the KLF with a power-law with 
a slope of j6=0.21 ± 0.03, (J3 =dlogN/dK). Since this result is 
based on a number of roughly 60 stars within 1 .5" of Sgr A*, 
the observed stars themselves cannot account for a significant 
extended mass component. It is therefore necessary to extrap- 
olate the KLF at the faint end. However, we assume that the 
observed stars trace the mass carrying population which can 
currently only be accessed theoretically by extrapolating the 
KLF towards its faint end. 

For the extrapolation of the KLF two different slopes are 
considered, /3 = 0.21 and a steeper one, /3 - 0.35, which fits 
the Bulge KLF (Alexander & Sternberg, 1999) as well as a 
model of an old (~ 10 Gyr) stellar cluster of solar metal- 
icity (Zoccali et al. 2003). In the following we concentrate 
on the region within Ri of the inner Plummer model compo- 
nent. The total number of stars in the inner cusp (R< 0.4" 



Table 3. Extinction free mass-to-light ratios, M/L(2fjm), for 
different stellar types in luminosity class V at the distance of 
the Galactic Centre. 
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Spectral types 


1.29 
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1.60 
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1.93 


17.005 


1.40 


F5 


2.73 


17.695 
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GO 


5.02 


18.665 


0.79 


KO 


5.36 


18.915 


0.67 


K5 


11.77 


20.065 


0.51 


MO 


9.23 


20.065 


0.4 


M2 


16.66 


20.915 


0.33 


M3 


38.49 


22.315 


0.21 


M5 


87.55 


23.815 


0.12 


M7 


166.43 


25.265 


0.06 


M8 



or 15.5 mpc) brighter than a given magnitude can be directly 
estimated from the (extrapolated) KLF. 

The KLF gives the number of stars per surface area per 
magnitude. In order to calculate the total number of stars and 
the total stellar mass present within the spherical volume en- 
closed by the core radius of R= 0.4" of our central Plummer 
model component, we de-projected the KLF. Tablel^lists the 
resulting numbers of stars, the total cluster mass, Mc/, the 
average stellar mass in the cluster, Ma,,e,., and the correspond- 
ing MjL {Ifim) for clusters given by KLFs with slopes of 
y6 - 0.21 and (3 - 0.35 and for different cut-off magnitudes, 
between K= 20 (M„„„ = IMq) and K^l'S, (M,„,„ = O.O6M0). 
Mc/ is calculated from the observed KLF using Aj^-li and the 
M/L (2//m) values listed in Tabled 

The numbers listed in Table|3]show that the mass of stars 
present within R - 0.4" deduced from a/? = 0.21 slope KLF 
cannot be higher than 800 M©, even after extrapolation to the 
faintest magnitudes (see Table |3 for a list of K-magnitudes 
andO.1051 corresponding MS-stars at the GC). Therefore, a 
cluster with a KLF ofyS - 0.21 cannot explainamass of3100 
Mq within R= 0.4", estimated from direct number density 
counts with an M/L{2fi m) = 2 Mq/Lq (see § 2). 

M/L (2yum) converges to a value 2 Mq/Lq foryS = 0.35, 
and a magnitude limit between 26 and 27, in this case we 
find the value of 3100 M© required by our mass distribu- 
tion modelling in §2. Lower magnitude limits would increase 
M/L (2/im) further, and the mass of the inner cusp calculated 
from our Plummer model would be underestimated. Stellar 
types corresponding to the required M/L (2yum) ratios are 
listed in Table|3] 

Table|3]also shows that for a main-sequence stellar popu- 
lation M/L(2/im) does not exceed 4.0 M0/L0, even if we con- 
sider very faint/low-mass stars {K - 28) and a steep faint-end 
KLF (/? - 0.35). This value corresponds to a cusp represent- 
ing only 0.2% of the total mass (see Table|5}. We conlude that 
the amount of extended mass allowed inside the core radius 
of the 'inner cusp' would be far too large to be explained by 
a stellar population of MS stars, and therefore would require 
another type of mass carriers. We show in § 5.3 that such a 
configuration is possible. 
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Table 4. For each 'observed' K-magnitude limit (not corrected for extinction) we give the numbers of stars A^, cluster mass 
Ma, average stellar masses Maver, and mass-to-luminosity ratio, M/L{2fim), within the central 0.4" for hypothetical Plummer 
model type clusters with KLFs extrapolated by two power-law slopes of [} - 0.21 and (3 - 0.35. Mc/ is the mass inside the 
core radius Ri of the innermost Plummer model component. The derived MjL values include an extinction correction for 
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Table 5. Mass -to-light ratios M/L(2jum), for the different ojf- 
set position fitting orbits obtained from our model calcula- 
tions. The errors on M/L are model dependent. From scaling 
to the enclosed mass estimate in Fig.|2lwe estimate an error 
of 5%. Columns 3, 4, and 5 give the cusp mass inside the or- 
bit of S2 for core radii of respectively 20.2mpc, 15.5mpc and 
13.2mpc. 
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0.104 
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5.2. Stability of a Cluster of Low-Mass Stars 

For a multi-mass stellar distribution, high mass stellar rem- 
nants (stellar black holes and/or neutron stars) are expected 
to migrate to the centre as a consequence of dynamical fric- 
tion. One would expect that, within a Hubble time, these com- 
pact objects show a higher concentration toward the centre 
than the lighter ones (Morris et al., 1993; Miralda-Escude 
& Gould, 2000), which should be transferred by that mech- 
anisms to orbits at greater distances from the centre of the 
cluster. This argues against the existence of a cluster of low- 
mass stars in the inner cusp. 

On the other hand, the possibility of a cluster of low-mass 
stars cannot be excluded and we are far from understanding 
the properties of the stars in the cusp. There is for example the 
unexplained presence of massive young stars, e.g. MS O/B- 
type stars close to the black hole (Genzel et al., 1997; Eckart, 
Ott & Genzel, 1999; Figer et al., 2000; Gezari et al., 2002; 
Ghez et al., 2003). These stars have not had enough time to 



achieve energy equipartition with the fainter older stellar pop- 
ulation. They are hence dynamically un-relaxed. Also, there 
are indications for a radial anisotropy of the stars in the cusp 
which might be un-relaxed (Schodel et al., 2003), in spite of 
the expected short relaxation time in this dense environment. 
Because of this general lack of theoretical understanding of 
the cluster near Sgr A*, we consider that a cluster of faint/low 
mass stars should not be ruled out entirely from the possible 
interpretations of the inner cusp. 

5.3. Is the Cusp Dominated by Dark and Massive 
Objects? 

In this section, we consider orbital fits with errors of < 1 cr that 
result in high cusp masses yet are in good agreement with 
the enclosed mass measurements. The evaluation presented 
in § 5.1 shows that such a heavy cusp is unlikely to consist of 
stars only. Here we study whether such a cusp could consist 
of heavier mass carriers like stellar black holes or neutron 
stars. 

5.3.1. Presence of stellar remnants in the center due to 
dynamical friction 

When the compact remnants of massive stars are themselves 
significantly more massive than the normal field stars in the 
Galactic Centre, as would be the case of black hole remnants, 
than they are susceptible to inward migration as a conse- 
quence of dynamical friction. The resulting mass segregation 
can lead to a pronounced concentration of compact objects 
within the central stellar core within a Hubble time (Mor- 
ris, 1993). A^-Body simulations of globular clusters showed 
that the combination of stellar evolution (production of stel- 
lar mass black holes, neutron stars, white dwarfs, and of bina- 
ries including such objects) and stellar dynamics will almost 
certainly lead to a strong increase of Mj L in the central parts 
of the nuclear star cluster Black holes and neutron stars sink 
to the centre and may coalesce (Giirkan et al. 2003). Such 
highly detailed self-consistent simulations of the dynamical 
episodes are, however, not yet possible on the scale of the 
Galaxy. A^-Body simulations (even with cutting-edge special 
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purposes computers like GRAPE-6 (Makino 2001)) cannot 
follow the evolution of a galactic nucleus over a Hubble time 
if relaxation is appreciable (Freitag & Benz, 2002). 

Under certain assumptions and for some range of initial 
parameters, less realistic but more efficient Monte Carlo nu- 
merical simulations of the evolution of the GC (Freitag & 
Benz, 2002) showed that the SBHs sink to the center on a 
short timescale of a few gigayears, settle into a centrally con- 
centrated distribution and dominate the stellar mass there. 
Similarly, Murphy Cohn & Durisen (1991, see their figure 
8b) showed that densities higher than 10^ M© pc"-' could re- 
side in regions as close as few mpc from the central SBH. 
Other studies by Morris (1993) and Miralda-Escude & Gould 
(2000) estimate that 10"* - 10^ SBHs, due to dynamical fric- 
tion, would have accumulated at distances of less than about 
1 pc from the center. These latter works can only account for 
about 10^-10"* Mq within separations of about 20 mpc from 
the centre which represents only 1% to 10% of our upper 
mass limit derived from the orbit of S2. However, the cen- 
tral density of the SBH cluster depends on various uncertain 
quantities: the SBH mass function, the stellar IMF and forma- 
tion rate, the remnant progenitor masses, and the dynamical 
age of the GC. Morris (1993) argued that within a wide range 
of assumptions about the IMF, and about the minimum stel- 
lar mass capable of producing a black hole remnant, the total 
mass of remnants concentrated into the inner few tenth of a 
parsec, would be about 0.4-5 10* Mq. If the black hole rem- 
nants were to achieve equipartition with the field stars, they 
would form an iimer core with a radius of a size as smaU as 
50 mpc. These would coalesce or form a quasi-stable clus- 
ter of stellar mass black holes. Even if a catastrophic merger 
of stellar remnants did occur at some point at the age of the 
galactic nucleus, the continuous influx of massive stellar rem- 
nants would ensure that a concentrated population of them is 
present within the steUar core. These findings apply, however, 
in the absence of an initial central black hole. To our knowl- 
edge, a similar study of a cluster with already a pre-existing 
super-massive black hole, as in the centre of our Galaxy, is 
not yet available. 

5.3.2. In situ formation of stellar remnants 

In the case of the Galactic Centre, complex dynamical 
episodes have taken place. A numerous variety of young 
early-type, bright and massive stars exist at distances of 
10 mpc - 400 mpc from the centre. There exist a dozen 
of bright 0/B main sequence stars within about 40 mpc of 
Sgr A*, these are fast moving S-stars (mass ~ 20 Mq) of 
which the star S2 is an example. There exists also ~ 30 more 
massive (30 Mq < mass < 100 Mq) very bright early-type 
stars, the so-caUed He stars exhibiting He/Hi emission lines. 
The existance of these two types of stars at these regions 
from the center is still not understood. These could have 
formed there or would have migrated there due to different 
in-spiraling processes (Genzel et al. 2003; Ghez et al. 2003; 
Gerhard et al. 2001, Krabbe et al. 1995). If such stars were 
always present in the center of the Galaxy, due to stellar evo- 
lution, their remnants would contribute very efficiently to the 



formation of a dense cluster also at the center. Here, we as- 
sume that these stars would end their Uves in the formation 
of neutron stars or stellar black holes with masses between 
1.5 Mq and 10 M© (on average 7 Mq), we also consider their 
lifetimes to be < 10^ yrs. If after a single Ufetime about 40 
of such stars form stellar remnants, throughout the age of the 
Galaxy (~ 10^" yrs), it is possible to account for about 10^ M© 
needed to explain our upper mass hmit of the 'inner cusp' . 

Considering the two above described scenarios - while 
these do not represent the complete history of the galactic 
nucleus - it is, however, a fair conclusion that strong to mod- 
erate dynamically caused M/L variations prevail at the GC. 

Our attempts in this work is not to explain the formation 
of a high density cluster of stellar remnants at the GC. We 
would rather like to analyse such a configuration if existant. 
In the following we investigate if such a hypotheticaUy high 
M/L configuration of steUar renmants can form a stable 
configuration. 



5.3.3. Stability of a cluster of stellar remnants 

Rauch & Tremaine (1996) studied the configuration of a cen- 
tral massive black hole plus an extended mass distribution M 
of radius R consisting of objects with mass m in terms of its 
non - resonant relaxation time f"^^. Under the assumption that 
M « Mbh, Rauch & Tremaine (1996) derive how f"^, de- 
pends on Mbh, M, and the orbital time scale tgrb at the outer 
edge of the cluster If (a) the stellar orbits have random ori- 
entations and moderate eccentricities, (b) the density of stars 
is approximately uniform within R, and (c) Mduster I Mbh 
1(1 - 10"^ then Rauch & Tremaine (1996) find that 



Mi 



n^NlnA 



(5) 



where IriK is the Coulomb logarithm ~ 13 in this case. 
This situation should - to first order - be applicable to the 
Black Hole/cusp scenario at the Galactic centre. Condition (a) 
is probably fulfilled with the possible exception that the stars 
in the cusp might have fairly high eccentricities (Schodel et 
al. 2003). Assuming a Plummer model as a cusp descrip- 
tion fulfills condition (b). We consider the case where a peri- 
centre-shift of the order of 1° is induced. As an example, we 
choose the case of a cusp mass of 15% with a total mass of 
4.1 X 10* Mq. This cusp mass is well below the upper hmit 
derived in section 4.1. In the case of a cusp mass of 15%, 
Ml Mbh is of the order of 10"^, if we restrict ourselves to the 
region within the core radius Ri. This is close to what is re- 
quired by condition (c). 

Here, Mbh = 3.5 X 10* and nisr is the average mass 
of the SteUar remnants. A value of rrisr = 5 M© is roughly 
consistent with a composite cluster made of neutron stars (m 
~ 1.5 Mq) and stellar black holes (5 Mq < m < 25 Mq). We 
then find that t^^^ ~ 10* X torb, if we assume that most of the 
mass inside the core radius is present in the form of stellar 
remnants of average mass rusr- Considering that the core 
radius of 15-20 mpc will define torb we find t"^^ to be about 
2 X 10^ years. 
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We can also investigate whether such a configuration is 
stable by estimating how many stellar black holes evaporate. 
Integrating a Maxwell distribution function for the velocities 
above the escape velocities gives the percentage of stars not 
bound to the BH. For the velocitiy dispersion, following 
Alexander et al. (2003), we can write: 

Vescape = ^Vcircular - ^(2(1^^ CT. (6) 

For a value of a = 2, about 0.03% of the present stellar 
black holes will be evaporated after each relaxation time. For 
a steeper cusp with a-3 we find that 0.006% of the stellar 
black holes will evaporate on that time scale. For f"^^^ of a few 
10^ years and a between 2 and 3 about 50% of the stellar 
black holes will have evaporated after 25 to 250 f"^, corre- 
sponding to about 5 X 10'^ to 5 X 10^ years. Given the fact that 
there will also be an influx of mass i.e. stellar black holes or 
neutron stars from outside the cusp one can consider such a 
configuration stable over a significant fraction of the Milky 
Way's age. 

Given equation (|5jl, for a fixed torb, the relaxation time 
decreases linearly with A^, i.e with the cusp mass. Thus, cusp 
masses well exceeding 2x10^ M© would not form a stable 
configuration compared to the age of the Milky Way. 

6. Summary 

From stellar number counts, we have a strong reason to as- 
sume that there is some extended mass around Sgr A*. There- 
fore, orbits of stars near Sgr A* should be analyzed in terms 
of non-Keplerian motion. Only in that way it is possible to 
find out more about the composition of the cusp. This paper 
presents the first approach to this problem. In this work we 
studied the case where the contribution of the central stellar 
cluster at the centre of the Milky Way could introduce peri- 
centre-shifts of the order of 1 degree per revolution on the 
orbit of S2. We assume that the dynamics of S2 is governed 
by both the super-massive black hole and the central cluster 
The steeply rising, cusp-like stellar density distribution in the 
inner few tens of milli-parsecs near Sgr A* can be modeled 
with a Plummer model distribution having a relatively small 
core radius of the order of the S2 semi-major axis. In a spher- 
ical potential, the stellar orbits will precess and the stars will 
follow rosetta shaped trajectories (Binney & Tremaine, 1987; 
Rubilar & Eckart, 2001; Fragile, & Mathews, 2000). We per- 
formed orbital fits to the S2 data for this class of orbits. We 
used a fourth order Hermite integrator, adequate for the re- 
quired precision. The runs were performed for a large range 
of initial parameters for the total mass, the position of the 
centre of mass, the fraction of the extended mass component, 
and the 6 classical orbital parameters. 

We found that the S2 trajectory can be equally well fit- 
ted by non - Keplerian orbits as by Keplerian orbits. For the 
Keplerian case, we could confirm the results of Schodel et al. 
(2003), Ghez et al. (2003), and Eisenhauer et al. (2003). Ba- 
sically independent of the amount of total or extended mass 
we find consistently a value for the central Black Hole mass 
of about 3.7 x lO^M©. This is well between but not coinci- 
dent with one of the two possible Black Hole mass estimates 



of either 4.75 x IO^Mq or 2.79 x 1O''M0 that were derived 
by Aschenbach et al. (2003) from an analysis of the quasi- 
periodicity observed in strong X-ray flares. 

A non - Keplerian orbit with a fit error of < 1 cr that 
shows a peri-centre-shift of the order of 1° implies a total 
mass of 4. 1 X 10* M© with a fraction of the extended mass be- 
ing / = 0. 1 5 . With this value of /, the amount of mass present 
in the inner cusp is ~75 times larger than the one expected 
form a direct mass estimation deduced from the stellar den- 
sity and the integrated mass plot ( see Fig.^and Fig.|2li while 
assuming a constant 2/im mass-to-light ratio of 2 Mq/Lq. 

If the inner cusp is mainly composed of main sequence 
stars, then a MjL (2/im) of 2 Mq/Lq can only be attained in 
the case of a faint-end KLF slope of yS = 0.35 and a magnitude 
limit 25 < A" < 26. The mass of such a cluster will be 3 100 M© 
within R - \5 mpc derived from our initial Plummer model 
in § 2. The existence of even fainter {K < 28), low-mass stars 
(0.06 Mq) in the inner cusp will not contribute more than 
~3 times the mass in the inner cusp; such a configuration is 
not likely to be stable for a long time. In addition it would 
be unclear how these stars would get into the inner cusp. It 
is unlikely that they formed there, however, they might have 
migrated into that region by loosing angular momentum by 
interactions with stars in the overall Galactic Centre stellar 
cluster with a core radius of 0.3-0.4pc. 

The example of a total mass of 4.1x10''Mq and an ex- 
tended mass component of 15%, as discussed above, implies 
an M/L(2jum) of ~ 145 + 10 Mq/Lq ; here objects with large 
MjL values are needed to explain the extended cusp mass. 
One plausible explanation could be the presence of a clus- 
ter of massive stellar remnants. By applying equation. |5] we 
find that an equilibrium configuration is possible under cer- 
tain conditions. 

In this study relativistic effects were not included for the 
orbit of the star S2. The peri-centre-shift induced by rela- 
tivistic effects is of the order of 9 arcmin compared to the 
40 arcmin expected from the Newtonian shift due to the pres- 
ence of an extended mass (see also Rubilar & Eckart, 2001, 
Fragile, & Mathews, 2000). It should be pointed out that the 
peri-centre-shift due to relativistic effects is in the opposite 
direction of that due to the extended mass and that relativis- 
tic effects could in a way mask the amount of extended mask 
around the black hole. Therefore it will be important for fu- 
ture work to investigate relativistic effects in more detail. 

After only a few additional observational epochs it will 
be possible to put even stronger constraints on the mass of the 
cusp, and consequently to know more about the mass carrying 
population in the immediate vicinity of the black hole at the 
centre of the Milky Way. 

As discussed in Rubilar & Eckart (2001) a complete de- 
termination of the cusp parameters can be expected from 
combinations of 3 or more orbits since then the number of 
degrees of freedom for the fits is increasing rapidly. Further 
studies with more general density distributions should be un- 
dertaken, because they may influence the precise value of the 
extended mass contribution. 
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Appendix A: Orbital parameters of three 1 cr 
fitting orbits 

In table A.l, we give the orbital parameters of 3 exemplary 
orbits giving fits within a 1 cr limit. The orbits corresponding 
to those fits are plotted in Fig.6. These are chosen for the case 
where the central black hole is present at the offset position 
((R.A - 0.09 mpc; Dec - -0.12 mpc) from the central nomi- 
nal radio position of Sgr A*). We also give the full range of 
parameters [Min;Max] which gives fits corresponding to a So- 
level. Here, we exclude solutions with M,o, > 5.4x10^ Mq. 
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Table Al. Orbital parameters of three exemplary 1 a fitting orbits. 



Total Mass M,„, (10 Mq) 


3.65 


4.10 


4.80 


Min 


Max 


3.28 


5.40 


roint Mass Mbh (10 Mg) 


i.o5 


3.69 


3.65 


3.28 


3.97 


Cusp Percentage(%) 


0.0 


10 


25 





35 


M"*^" : Cusp mass inside !s2 s orbit(10 Mq) 


0.001 


0.049 


0.133 





0.26 


e : bccentricity 


0.88 


0.882 


0.881 


0.868 


0.897 


(' : Inclination(°) 


-47.89 


-47.88 


-48.12 


-46.86 


-50.19 


u '. Semi-major axis(mpc) 


H.OJ 






A "nO 


A 70 


T : Time of peri-centre passage 


2002.327 


2002.328 


2002.333 


2002.327 


2002.337 


CI : Angle of line of nodes(°) 


45.3 


45.38 


44.59 


42.17 


46.1 


a> : Angle to peri-centre(°) 


145.1 


144.84 


145.6 


144.1 


147.7 


P : Period (Years) 


15.55 


15.43 


15.06 


14.61 


16.67 


Peri-Ccntrc distance 


0.557 


0.552 


0.55 


0.511 


0.6 


Aor : Peri-centre-shift per revolution 


OO^OO'OO" 


00°40'48" 


01°52'43" 


OO-OO'OO" 


03''07'12" 
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